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The low-momentum interaction Vj ow -k derived from realistic models of the nucleon-nucleon inter- 
action is presented in a separable form. This separable force is supported by a contact interaction in 
order to achieve the saturation properties of symmetric nuclear matter. Bulk properties of nuclear 
matter and finite nuclei are investigated for the separable form of V[ ow _k and two different param- 
eterizations of the contact term. The accuracy of the separable force in Hartree-Fock calculations 
with respect to the original interaction Vi ow -k is discussed. For a cutoff parameter A of 2 fm _1 a 
representation by a rank 2 separable force yields a sufficient accuracy, while higher ranks are re- 
quired for larger cut-off parameters. The resulting separable force is parameterized in a simple way 
to allow for an easy application in other nuclear structure calculations. 
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I. INTRODUCTION 

The evaluation of bulk properties of finite nuclei and nuclear matter starting from realistic models of nucleon- 
nucleon (NN) interaction is a major challenge in modern nuclear physics. Since the exact form of the interaction 
resulted from the underlying theory of the strong interaction remains unknown, one usually has deal with realistic 
models developed so as to fit experimental data for free nucleon-nucleon scattering up to the threshold for pion 
production and properties of the deuteron [J-Q- It was done by obtaining a best fit for a large number of adjustable 
parameters using several thousands experimental points so that there exists several quite different potential models 
commonly used. A general feature of such realistic interactions is strong short-range and tensor components, which 
- 1 . cannot be handled within the standard perturbation theory. There have been suggested different approaches in order 
to overcome this problem: Bethe-Brueckner-Goldstone expansion k5|, correlated basis functions [6j, quantum Monte 
Carlo Q, self-consistent Green's function theory (see, e.g., Ref. Q). These methods were successfully applied to 
describe bulk properties of nuclear matter pairing gap of nucleons [Ioj|, weak response [ll| and shear viscosity 
fNj ', of nuclear matter [l2j. However these approaches remain very complex to be applied directly to a description of 
finite nuclei, as well as inhomogeneous nuclear matter, also known as pasta phase, which exists in the inner crust of 
neutron stars. Alternatively, they have been combined either to phenomenological approaches through a local density 
£f) i approximation [l3j . or as an input for a density functional approach [Til ]. In these approaches, adjustable parameters 
' need however to be determined. 

Besides the realistic interactions, various phenomenological models have been developed, such as the Skyrme inter- 
action [l5| . and adjusted to describe the experimental data for the ground states of finite nuclei and the empirical 
saturation point of symmetric nuclear matter. A simple parameterization of such phenomenolo gica l forces through 
the local single-particle densities allows a simple solution of the Hartree-Fock (HF) equations [la]. Finally, these 
models have been successfully used for predictions of equations of state (EoS) of nuclear matter and description of 
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■ pasta phase within the Wigner-Seitz (WS) cell approximation 17[. In neutron stars, these models are extrapolated 
far from the condition where it has been adjusted and might in some cases become unstable (l8j . The instabilities of 
these models could however be corrected such as it reproduces the features of a G- matrix in nuclear matter [l9j . 

An alternative method, which is based on realistic NN interactions and allows us to perform Hartree-Fock calcu- 
lations similarly to the phenomenological forces is the low- momentum interaction Vow-k- The basic idea of Viow-k 
is to separate the predictions for correlations at low momenta, which are constrained by the NN scattering matrix 
below the pion threshold, from the high-momentum components, which may strongly depend on the underlying model 
of realistic NN interaction. By introducing a cutoff A in momentum space, one separates the Hilbert space into a 
low-momentum and a high-momentum part. The renormalization technique (see, e.g., [20l - [24l |) determines an effective 
Hamiltonian, which must be diagonalized within the model space (below the cutoff). With the cutoff in the range 
of A = 2 fm _1 Vow-k becomes model independent, and reproduces the deuteron binding energy, low-energy phase 
shifts, and half-on-shell T matrix with the same accuracy as the initial realistic interaction. This model independence 
demonstrates that the low-momentum physics does not depend on details of the high-momentum dynamics. 

In spite of its obvious advantages V ow -k potential still remains a quite complicated object. On the one hand, it is 
nonlocal and therefore is represented as a matrix element in momentum space for each partial wave channel. This 
nonlocality increases the computational time in Hartree-Fock itcractions, and prevents the use of Viow-k if the number 
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of nuclcons is too large, such as in the Wigncr-Scitz cells present in the crust of neutron stars for instance [17|. On the 
other hand, the renormalization technique used to produce Vjow-k seems not to be trivial. The resulting interaction 
is given as a matrix table which is not an easy-to-use form and prevents this potential to be popular. A possible 
way out is to find a separable representation of Mow-kj since it significantly simplifies many-body calculation |25l. [26|. 
Moreover recent calculations of triton binding energies demonstrate the Vi ow -k can be very good approximated by a 
low-rank separable force for low values of the cutoff A [23]. We investigate the separability of Vow-k by using the 
diagonalization of the matrix in momentum space for each partial wave channel. It allows us to find a low rank 
separable form of Vow-k, which can be used in HF calculations of nuclear matter as well as finite nuclei. 

The Viow-k Hartree-Fock calculations demonstrate a monotonic increase of the binding energy of symmetric nuclear 
matter as a function of the nucleon density, thus it cannot reproduce the empirical saturation point (2a . [30| . Therefore 
we supplement Vow-k by a simple density-dependent contact term, which accounts for a three-body correlations. This 
contact term is adjusted to reproduce the saturation property of symmetric matter. 

The paper is organized as follows. In the next Section [TT| we discuss the model space technique used to produce 
Vow-k and outline the procedure to determine the separable representation. In the last Section IIIII we sum up all 
results and suggest a simple fit for separable representation of V\ w-k as well as two different parameterizations of the 
contact term, adjusted for the fitted potential. 

II. MODEL OF THE NN INTERACTION 

The main idea of Vi ow -k interaction is to disentangle the low-momentum or long-range part of a realistic NN 
interaction, which is fairly well described in terms of meson-exchange, from the high-momentum or short-range part 
where quark degrees of freedom are getting important. In other words, one defines a model space, which accounts for 
the low-momentum degrees of freedom and renormalizes the effective Hamiltonian for this low-momentum regime in 
order to account for the effects of the high-momentum components, which are integrated out. 

In practice the Vjow-k interaction can be derived either using model space methods (such as Lee-Suzuki [2l| or 
Okubo [22| ) or through an renormalization group treatment [2(j ■ Both approaches are essentially equivalent and lead 
to the same energy- independent potential [23(. In the following we will use the model space technique to disentangle 
these parts based on the unitary model operator approach (UMOA). This approach has frequently been described in 
the literature [HI, [28|, Therefore we will restrict the presentation only to basic equations, which will define the 
nomenclature. 

To determine the model space, the low-momentum subspacc of Hilbcrt space, one defines a projection operator P, 
which projects onto this model space. The complement of the subspacc will be defined by the projection operator 
Q, in such a way that the whole space is covered by these two operators. Thus they satisfy the following relations 
P + Q = 1, P 2 = P, Q 2 = Q, PQ = = QP. The unitary model operator approach defines a unitary transformation 
U in such a way that the transformed Hamiltonian does not couple the P and Q space, i.e., 

QU~ 1 HUP = 0. (1) 

Now the effective two-body interaction of Hermitian type can be determined in terms of unitary transformation U as 

V eS = V M = U- 1 (h -v 12 )U-h , (2) 

where #12 stands for the bare NN interaction. The operator ho denotes the one-body part of the two-body system 
and contains the kinetic energy of the interacting particles. It is im por tant to notice that in any case ho commutes 
with the projection operators P and Q. As it was shown by Suzuki [24j the operator U is expressed as 

U ={1 + lo-lo v ){1+CjCj^ +^Cj)~ 1/2 1 (3) 

where an operator ui fulfills relations u> = QCjP and ui 2 = Cj^ 2 = 0. To evaluate the matrix elements of this operator 
w one should first solve the two-body eigenvalue equation 

{ho+vi2)\$k)=E k \<Z> k ), (4) 

From eigenstates \$k) we determine those eigenstates |$ p ), which have the largest overlap with the P space. After 
the respective matrix elements of w and later U may be defined in terms of P{Q) eigenstates. This matrix element of 
U can then be used to determine the matrix elements of the effective interaction V c g in P space (for the details see 
[HI [HI)- I n this way, one obtains the effective Hamiltonian H c r = ho + V c ff. Diagonalising it in the low-momentum 
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FIG. 1: (Color online) Top: Nonzero eigenvalues ch of 1 Sq channel. Bottom: Nonzero eigenvalues aj of 3 Si- 3 D\ 

channel. 

model space (P space) , one obtains eigenvalues which are identical to the diagonalization of the original Hamiltonian 
ho + V in the complete space. Moreover, the solution of the Lippmann-Schwinger equation for NN scattering phase 
shifts using Vi ow -k with a cutoff A yields the same phase shifts as obtained from original interaction V12 without a 
cutoff. If the underlying interaction is a realistic interaction, fitted to reproduce the experimental phase shifts below 
A, these phase shifts will be also reproduced by Viow-k- 

If the cutoff A is chosen around A = 2 fm _1 the resulting T4 ow -k is found to be essentially model independent, i.e., 
is independent on the underlying realistic interaction vi2- In this sense I4ow-k is unique and, as it reproduces the NN 
scattering phase shifts it can also be regarded as a realistic interaction as, e.g., the CD-Bonn [l[ or Argonne V18 Q 
potentials. 

Originally Vj w-k is nonlocal and defined in terms of matrix elements in a basis of NN states labeled by relative 
momentum for pairs of nucleons. Thus for each partial wave channel there exists a matrix, which represents Vjow-k {k, k') 
on a mesh of N discretized relative momenta k and k' in the range < k, k' < A. Since this matrix is real and 
symmetric with respect to k, k' one can diagonalize it, so that, it can be written as a sum of N real eigenvalues 
multiplied with the respective eigenvectors 

TV 

V M {k,k') = Y,aif*{k)fi{k'), (5) 
t=i 

where N is the number of mesh-points and the dimension of Viow-k matrix. The eigenvectors fi{k) satisfy the 
orthogonality relation 

- [ dkk 2 f i (k)f j (k) = 5 ij . (6) 
n Jo 

In the following we will omit the symbol of complex conjugation because all eigenvectors are real. The last equality ([5]) 
is nothing else but the general definition of a separable potential of the rank N. If the rank of the separable potential 
equals the dimension of the matrix Vi ow _k(fc, k') the whole information is exactly restored from the eigenvalues and 
eigenvectors /j. As we will see later, some of eigenvalues can be zero or negligibly small so that one can reduce 
the rank of separable interaction taking into account only the n eigenvalues with largest absolute values. It leads to 
a new approximated separable interaction V^_ k (fc, k') 

n 

V low . k (k, k') ~ V™_ k (k, k') = ]T (HMmW), (n < N). (7) 

i=l 

The low-rank separable representation of NN interaction leads to significant simplifications in many-body calculations. 
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FIG. 2: (Color online) Squared deviation of the separable V£j_ k (fc, k') from the original Vi ow _k(fc, in 1 Sq channel 

for different values of the cutoff parameter A. 



The effective interaction Viow-k as well as its separable form is nonlocal and defined in terms of matrix elements in 
momentum space. It implies that the HF calculations has to be performed in a Hilbert space using an appropriate 
basis | a), | /?},.. . The HF Hamiltonian is then expressed in terms of the matrix elements between these basis states 
(a | H^f I P) and the HF single-particle (s.-p.) states | ^> n ) are expressed through the expansion coefficients in the 
basis 

\^ n ) = J2\a)(a\^ n )=J2^\a). (8) 

a a 

The part of the HF Hamiltonian originating from Vjow-k can be expressed in terms of two-body matrix elements by 

{a | #iow-k|/3) = ^<a 7 |Mow-k|/?<5}p 7 a, (9) 

7,6 

where p y s is the single-particle density matrix. In order to investigate the bulk properties of finite nuclei we perform 
HF calculations within the spherical Wigncr-Scitz cell assuming a plane wave single-particle basis [3lL l32j. 



III. RESULTS AND DISCUSSION 



In the following we discuss results for symmetric nuclear matter as well as finite nuclei obtained from HF calculations. 
These calculations arc performed in the model space, which is defined by a cutoff parameter A in the two-body 
scattering equation, employing the corresponding low- momentum interaction Vjow-k, which is derived from the CD- 
Bonn [l[ interaction using the technique described in Sec|Hj The NN interaction has been restricted to partial waves 
with total angular momentum J less equal 6. 

We start our discussion with the comparison of the eigenvalues ai obtained from diagonalization of 20 x 20 matrix 
of Mow-k(fc,fc') in 1 5o channel. The resulted nonzero eigenvalues are shown on the top of Fig Q] for different values 
of A. As it was discussed above, Viow-k interaction becomes model independent at A = 2 fm -1 . At this value of the 
cutoff parameter A the diagonalization procedure yields only 11 nonzero eigenvalues, other words, Mow-k interaction 
in 1 Sq channel is a separable interaction of the 11th rank or, following ([5]), one can write 

^(MO = Mow-k(M')- (10) 

The nonzero eigenvalues are essentially independent on N, the dimension of the matrix representing Vj ow -k- Going 
further one can notice, that many of the nonzero eigenvalues are nevertheless very small, and only some of them, e.g., 
at i = 1,2,20 carry the main part of the information about the interaction model. This gives rise to a substantial 
lowering of the rank of separable potential, as it was shown in Eq.©. With the increase of the cutoff A the absolute 

values of the eigenvalues increase as well and as a consequence the rank n of the separable form Vj[^J_ k defined in (J7J has 
to be increased to achieve a reasonable accuracy. Increasing A more information about the short-range components 
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FIG. 3: (Color online) Squared deviation of the separable V^u(k,k') from the original V\ ow -k(k, k') in 3 5i- 3 Z?i 

channel for different values of the cutoff parameter A. 

of the underlying bare interaction is included, which requires a larger rank in the separable representation of the 
interaction. 

In case of the coupled channels, like 3 Si- 3 Di channel, the dimension N of the Vjow-k matrix is twice as large if one 
keeps the number of mesh-points in each channel the same as for the uncoupled partial waves. It turns out that also 
the number of nonzero eigenvalues increases as shown in the lower panel of Fig[TJ It is obvious that the rank of the 
separable potential should be higher than for 1 S'o channel. It is a general feature that coupled channels require higher 
rank separable interaction than the uncoupled ones |33j ] . Also one observes pairs of positive and negative eigenvalues 
which have about the same absolute value. This picture remains for higher values of A. As we will see later, this 
symmetry between positive and negative eigenvalues will play a crucial role in the convergence of the separable form 

^Low-k t° the initial Mow-k with increase of rank. 

In order to determine a minimal rank for a reliable separable approximation in each channel we calculate the square 
deviation 77 of the separable form V^J_ k from the original potential Vjow-k for each rank n 



In Fig (21 the deviation for 1 Sq channel at different values of the cutoff A is shown. At A = 2 fm _1 one observes a fast 
convergence to zero deviation already at the rank n = 2. The growth of the cutoff monotonically increases the rank 
of the separable potential. At A = 3 fm _1 one may expect a good accuracy starting from n = 5. 

The deviation 77 for 3 Si- 3 -Di channel is displayed in Fig|3l First, at low n the absolute value of the deviation is 
one order of magnitude higher than for uncoupled 1 5o channel. Increasing the rank one observes a non-monotonic, 
oscillating decrease of 77, specially for high A. As we have seen, the diagonalization of the channel 3 5i- 3 Z?i yields 
both positive and negative eigenvalues, which are symmetrically distributed over i. So that they form "pairs" with 
very similar absolute values. Assuming the odd rank we take into account either uncompensated positive or negative 
eigenvalue. This eigenvalue will be compensated in the next (even) rank, and the accuracy will be significantly 
improved. 

The deviation 77 for various other channels at A = 2 fm _1 is shown in FigsJH [5J In the following we choose the 

second rank approximation for the uncoupled channels (n = 2) and the third rank for the coupled one (n = 3). Below, 

[2 31 

the respective separable version of Mow-k will be referred to as V[o W _ k - 

Now let us turn to the binding energy of symmetric nuclear matter, which is displayed in Fig. [5] The HF calculations 
using VJo W _ k (dashed line) yields essentially the same result as the one employing the original Vjow-k interaction (solid 
curve). The deviation does not exceed 1% at the saturation density po and 1.7% at the density 2po. We also compared 
the binding energy of pure neutron matter for both potentials and found that the discrepancy is less than 1% for the 
same range of densities. 

However, neither of the 2 calculations yields a saturation point, i.e. a minimum in the energy versus density plot, 
as it has been observed before [3(3, IH[ • This absence of the saturation is one of the main problems in calculations of 
nuclear matter employing Vjow-k- It cannot be cured by the inclusion of correlations beyond the HF approximation, 
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FIG. 4: (Color online) Squared deviation of the separable VjJ^l k (fc, k') from the original Viow-k(fc, k') for various 

uncoupled channels at A = 2 fm _1 . 
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FIG. 5: (Color online) Squared deviation of the separable V^_ k (k, k') from the original Vi ow -k(fc, k') for various 

coupled channels at A = 2 fin. -1 . 

e.g., by means of the BHF approximation [28j. Recent relativistic calculations by van Dalen and Miither demonstrate 
that saturation can be achieved within the 14 ow -k approach by inclusion of relativistic effects in dressing the Dirac 
spinors which are used to evaluate the underlying realistic interaction [35[ . 

[2 31 

All the results obtained so far indicate that ^ is an accurate low-rank separable representation of Vi ow -k inter- 
action. However, in order to make it accessible to other users, it should be parameterized in a simple form. Here we 
suggest the fitting function for all fi(k) in all channels 

fi(k) = a t + (ft exp fyfc 5 *) + Mi) sin(A:cri + A,), (12) 

which contains 7 parameters for each partial wave channel and each fi(k). In the Table HVl we summarized all 
parameters of the separable fitted form for uncoupled channels, while all parameters for the coupled channels are 
shown in the Table [V] By using the values from both tables one can reproduce the fitted version of V^ 3 ^ for a given 

[2 31 

partial wave channel. In the following we will identify the respective separable fitted potential as Vg t : . 

[2 31 

In order to check the accuracy of our fit we perform HF calculations of nuclear matter employing Vg, ' . The 
respective binding energy as a function of the density of symmetric nuclear matter are displayed on Fig IB] by a dashed- 
dotted line. One observes that at up to saturation density pq ~ 0.16 fm~ 3 the fitted potential V^' 3 ' reproduces the 
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FIG. 6: (Color online) Energy per nucleon of symmetric nuclear matter as a function of the density. Results of 

\i 31 

Mow-k interaction (solid line) compared with the separable form KLLi (dashed line) and the respective fitted form 
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(dashed-dotted line). 



Interaction 


t [McV fm 3 ] 


t 3 [MeV fm a+to ] 


x 3 


CT 


-584.1 


8330.7 


-0.5 


CT1 


-548.0 


7890.13 


-0.5 


CT2 


-565.467 


8180.0 


-0.5 



TABLE I: Parameters of the contact interaction defined in Eq. (|13p . The set CT was produced for V\ 
CT2 supply!/' 2 ' 31 . 



while CT1 and 



rn Ql 

results of VjLLu (red dashed), while at higher densities it becomes slightly less bound and lies closer to the original 

12 31 

Mow-k (solid). Thus the deviation of the fitted separable potential Vg t ' from Mow-k does not exceed 1% of binding 
energy. Not going into details we mention that the deviation rises mainly from 3 Si- 3 Di and ^P^-^F^ coupled channels. 

As we have already seen from FigjG] Vj w-k interaction as well as its separable form V^' 3 ' does not describe the 
empirical saturation point. In order to achieve the saturation in nuclear matter one has to add three-body interaction 
terms or a density-dependent two-nuclcon interaction. Therefore we support the low-momentum interaction by a 
simple contact interaction, which have been chosen following the notation of the Skyrme interaction [l5l [l6| 



Au = Av + Az/ 3 , 



(13) 



with 



A^ = fo [(2 + x )p 2 - (2x + l){pl + p 2 p )] 



(14) 



and 



A^ 3 



24 



t 3 p a [(2 + x 3 )p 2 ~(2x 3 + l)(pl+p 2 p )] 



(15) 



where p p and p„ are the local densities of nucleons while the total matter density is denoted as p = p p + p n . The 
values of a and xq were fixed at a = 0.5, xq = 0.0, while to, t 3 , x 3 were fitted in such a way that HF calculations using 
^ow-k or Vg t ' plus the contact term (|1 3|) reproduces both the empirical saturation point of the symmetric nuclear 



matter and the symmetry energy at saturation density. Following [32[ the contact interaction produced for Viow-k 

[2 31 

will be labeled by CT, and the respective interaction model Vi ow -k+CT. For the fitted potential Vg t ' we suggest two 
possible paramctcrizations: CT1 and CT2. Their parameters and properties of nuclear matter are shown in Tables U 
and|TTl respectively. 
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Density p [fin 3 ] 

FIG. 7: (Color online) Energy per nuclcon of symmetric nuclear matter as a function of the density. Results of 
Mow-k + CT interaction (solid line) compared with the fitted separable form V^' 3 ' + CT1 (dashed line) and Vg? + 

CT2 (dashed-dotted line). 



Interaction 


po [fm- J ] 


E/A(p ) [MeV] 


K [MeV] 


Vlow-k+CT 


0.16 


-16.0 


258 


V£' 3] +CT1 


0.16 


-16.1 


241.9 


Vg l t 2 ' 31 +CT2 


0.156 


-16.0 


240.5 



TABLE II: Bulk properties of symmetric nuclear matter derived from Viow-k and its separable representation. They are 
supplemented by the respective contact interaction. 

The interaction V^' 3] + CT1 gives the binding energy per nucleon of symmetric nuclear matter E/A=-16.1 MeV at 

the density po = 0.16 fm~ 3 . The HF calculations of nuclear matter (see FigfT]) for V^' 3 ' + CT1 give results (dashed 
line) very similar to the non-separable initial interaction Vjow-k + CT (solid line). However, in calculation of finite 
nuclei we observe a deviation of about 0.12MeV in the binding energy of light nuclei, like 16 (see Table |TTT|. The 
picture can be improved if we assume, that the saturation density is not defined exactly and allow for a small deviation. 
Along this line the second parameterization CT2 was produced. The interaction Vg ' + CT2 gives E/A=-16.0 MeV 
at the density po = 0.156 fm~ 3 . This corresponds to a small shift of the saturation point with respect to the initial 
Vi ow _k interaction (see FigJ?]). It allows us to improve the accuracy in the binding energies of finite nuclei: one can 
notice that the contact interaction CT2 leads to a better description than CT1. However, comparing the rms charge 
radii of nuclei in Table iHfl we see that due to the shift in saturation density the interaction V^' 3 ^ + CT2 yields larger 

radii than the interaction Vg t ' + CT1. 

For all models considered here the compressibility modulus at saturation density is in the range 240.5 < K < 258 
MeV. This means that the respective equations of state displayed in FigJ7]are rather soft, at least at densities up to 
about two times saturation density. Such a prediction of a soft EoS is in agreement with data extracted from heavy 
ion reactions. For example, heavy ion data for transverse flow [36] or from kaon production [U support the picture 
of a soft EoS in symmetric nuclear matter. This value for the compressibility modulus is also in agreement with that 
of the Skyrme interaction which reproduce correctly the breathing mode in nuclei (giant isoscalar resonance) [42j ■ 

IV. SUMMARY AND CONCLUSION 

In the last decade it has become popular to perform nuclear structure calculations using the low-momentum NN 
interaction Viow-k (sec for instance the recent Ref. [HI, H3|). This interaction is constructed from a realistic NN 
interaction by introducing a cutoff A in the relative momenta of the interacting nuclcons. We used a model space 
technique on the base of the unitary model operator approach to separate the low-momentum and high-momentum 
parts of the initial CD-Bonn interaction. The cutoff parameter A was fixed at A = 2 fm _1 so that a Vi ow -k is obtained, 
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Interaction 


u 


Oa 


Oa 


Oa 


2UiWm 




E/A [MeVJ 










now-k+ ( -' y 


7 G1 


Q ^7 


-0.4Z 


7 7^ 


7 7C. 


Vjif' 31 +CT1 


-7.79 


-8.56 


-8.35 


-8.78 


-7.76 




-7.84 


-8.58 


-8.37 


-8.79 


-7.76 


Experiment 


-7.98 


-8.55 


-8.67 




-7.87 




r c h [fm] 










Vlow-k+CT 


2.79 


3.50 


3.54 


3.68 


5.51 


l4 2 ' 31 +CTl 


2.81 


3.51 


3.55 


3.68 


5.52 


yV2,3] +CT2 


2.82 


3.53 


3.58 


3.71 


5.56 


Experiment 


2.74 


3.48 


3.47 




5.50 



TABLE III: The binding energy per nucleon and rms charge radii of finite nuclei. Experimental data taken from Refs. [38|-[4ll]. 

which is essentially independent on the underlying bare NN interaction. 

The resulted Mow-k interaction is nonlocal and defined in terms of matrix elements in momentum space for each 
partial wave channel. This allows us to use a diagonalization method in order to express the matrix elements in a 
separable form. We investigate the separability in different channels with increase of the cutoff A. It was found that 
at A = 2 fm -1 the low-momentum interaction can accurately be approximated by a low-rank separable interaction. 
This separable interaction is parameterized to make it accessible for other nuclear structure calculations. 

A density dependent contact interaction is added to reproduce the saturation property of infinite nuclear matter. 
HF calculations using this interaction model also reproduce the bulk properties of finite nuclei with good accuracy. 
We demonstrate that this new separable representation of the Vio W -k interaction, reproduces the results derived from 
the original Vi ow -k with a high accuracy. 

This work has been supported by the European Graduate School "Hadrons in Vacuum, in Nuclei and Stars" 
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TABLE IV: Parameters of V£ 1 for uncoupled channels. See Eq. (fT2|l . 
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TABLE V: Parameters of v£' 3] for coupled channels. See Ea.(fT2|>. 
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